arXiv:nucl-th/0503050vl 17 Mar 2005 


Light Nuclei within Nuclear Matter 


M. Beyer S.A. Sofianos'^ N. Furutachi*^ S. Oryu^^ 

Institute of Physics, University of Rostock, 18051 Rostock, Germany 
Department of Physics, University of South Africa, Pretoria 0003, South Africa 

^Department of Physics, Tokyo University of Science, Noda, Chiba 278-8510, 

Japan 


Abstract 

We investigate the properties of ^He, ^He, ®He, ^Li and nuclei in nuclear matter 
of finite temperature and density. A Dyson expansion of the many-body Green 
function leads to few-body equations that are solved using the Integro-Differential 
Equation Approach (IDEA) and the Antisymmetrized Molecular Dynamics (AMD) 
methods. The use of the latter method allows us to trace the individual movement 
of the wave packet for each nucleon and the formation and disintegration of quasi¬ 
nuclei in a changing thermodynamical nuclear matter environment. 
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Nuclear fragments detected in a heavy ion collision at intermediate energies 
play a central role to gain information on the properties and the dynamics 
of nuclear matter under extreme conditions. One particular interesting exam¬ 
ple is the liquid gas phase transition usually related to multi-fragmentation 
into light nuclei [1,2]. For more and other experiments see e.g. [3,4,5,6,7,8]. In 
practical calculations of heavy ion collisions [9,10,11,12,13,14,15,16,17,18,19], 
formation of clusters is described by kinematical (or geometrical) constraints 
that allow certain nucleons to combine to larger clusters or larger clusters 
to disintegrate. In some simulations these constraints are supplemented by a 
probability to occur that is either related to some £t parameter or in more 
microscopic approaches to the nucleon nucleon cross section. Depending on 
the specific approach, e.g. statistical multi-fragmentation models [9,10], mi¬ 
croscopic transport models such as Boltzmann type [11,12,13,14,15], quantum 
molecular dynamics [16,17,18], or recent more intermediate type of event gen¬ 
erators [19], the algorithms that deal with fragmentation are rather educated 

^ Gorresponding Author: michael.beyer@uni-rostock.de, tel. -|-49 (381) 498 6773, 
fax. +49 (381) 498 6772 


Preprint submitted to Elsevier Science 


9 February 2008 





and powerful. Despite the decisive role of fragmentation, the exact conditions 
under which fragments are formed and which properties they have in such 
an excited nuclear many-body system are little known. In addition, even for 
a particular heavy ion collision at a given energy the formation of fragments 
could be at different stages under quite different conditions. 

The parameters that are introduced into the simulation to accommodate frag¬ 
mentation might be calculated microscopically. To motivate the objective of 
our present work we remind the reader how the conditions of deuteron forma¬ 
tion are achieved in a microscopic statistical description of heavy ion collisions 
such as the Boltzmann-Uehling-Uhlenbeck (BUU) approach. To this end an 
overlap integral is introduced that cuts momentum components and hence 
geometrically describes the effect of Pauli blocking. To be more specihc, for¬ 
mation of deuterons is only allowed if [11] 

J d^Qf (q + < ^cut (1) 


holds, where Pc.m. denotes the c.m. momentum of the cluster and q the relative 
momentum of the nucleons inside the wave function (j){q) of the deuteron. The 
Pauli blocking is achieved by the nucleon momentum distribution function 
/. The cut-off Pcut is a parameter that needs to be hxed. Originally it has 
been treated as a £t parameter with a typical value of ~ 0.2. It has 
been shown that it is possible to connect this parameter to a microscopic 
calculation [14]. Using in-medium deuteron wave functions leads to Fcut — 0.15 
not quite the value used to fit to experiments. It could be that the higher value 
of simulates also evaporation from larger clusters that is not present in the 
microscopic approach used in [14]. The physical interpretation of (1) is due to 
the Mott transition that occurs when the binding energy of the cluster vanishes 
because of medium effects. For an early paper on a microscopic calculation of 
the deuteron’s Mott transition see [20,21]. In the present work we proceed 
further along the pathway of Ref. [14] investigating the Mott transition for 
larger clusters up to using rigorous few-body equations that include the 
dominant medium effect, i.e. self energy corrections and Pauli blocking. This 
also extends previous investigations on the behavior of the nucleon deuteron 
reaction [14,22,23,24,25], the ^He/^H [26] and the a-particle [27,28] along this 
lines. These investigations were based on the Alt-Grassberger-Sandhas (AGS) 
equations [29] which were modified to include the effects of the medium. 

Because of the Pauli blocking the effective interaction between nucleons be¬ 
comes energy [E) dependent. To be specific we choose a temperature of T = 10 
MeV that is characteristic for heavy ion collisions and also relevant for neu¬ 
tron stars. In the present work we first endeavor to transform the P-dependent 
nucleon-nucleon interaction to an P-independent one using the Marchenko in¬ 
verse scattering method [30]. This paves the way to go beyond three- and 
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four-nucleon systems described by integral equations in momentum space by 
using configuration space formalisms. We employ here two such formalisms, 
namely the IDEA [31,32] and the AMD [33,34] methods. The former method 
is based on the assumption that the dominant correlations in the medium 
are still of two-body nature and thus one can obtain Faddeev-type equations 
for the corresponding two-body amplitudes. The AMD method, on the other 
hand, enable us to study changes in the density of the clusters embedded in 
the medium and obtain expectation values (r) for the center of the wave pack¬ 
ets for each nucleon within the cluster. In our studies we consider, the ^He, 
"‘He, ®He, ^Li and the nuclei. 

The Dyson expansion of the many-body Green fnnction is used to obtain 
thermodynamic cluster Green functions for the A-nucleon system embedded 
in nuclear matter. The cluster Green functions are evaluated for an uncor¬ 
related medium. Hence, they can be formulated as resolvents and respective 
resolvent equation can be derived for A quasi-nucleons that can be tackled 
with few-body techniqnes, however, modihed effective potentials and self en¬ 
ergies (masses). At a two body level this leads to the well known t-matrix 
eqnation [35], 

T 2 (z) = 14 + 14(1 - /l - f 2 )Ro{z)T 2 {z). (2) 

with Ro{z) = {z — iJo)~^-The corresponding Schrodinger-type eqnation is 

(Ho - ;^)T(12) + ^(1 - /i - /2)14(12, T2')T(T2') = 0. (3) 

1 ' 2 ' 


The fi is the Fermi fnnction given by 




1 


( 4 ) 


with f3 = l/ksT. We have used an effective mass approximation, i.e. the 
quasi-particle energy e is given by 


2m 




2m, 


eflf 


+ ^0 


HF 


( 5 ) 


The constant shift Aq^ can be absorbed in a redehnition of the chemical 
potential peff = k' ~ Aq^. For simplicity we take the effective masses used in 
previous calculations [14]. Assnming that the two-body system is at rest in 
the medinm, then the kinetic energy term (3) is written as 


kj ^ ^ ^ 1 ^ 

2mefr 2mefr 2mefr 2 


( 6 ) 
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while the effective interaction that appears in (2) and (3) is 


V{E,r) = {l-2f{E))V2{r), 


( 7 ) 


where E = z — -Econt is the two-body center of mass energy, i^cont denotes the 
continuum edge as explained, e.g., in Ref. [20,21]. 

To transform this potential to an equivalent i?-independent one, we employed 
the Marchenko inverse scattering method [30] which does not require to choose 
a priori the shape and range of the potential. The Schrodinger equation is 
transformed to an integral one for the non-local function K (r, r '), 

OO 

Ki{r, r') + Ee{r, r') + J Ke{r, s)Ee{s, r')ds = 0 , (8) 


where the kernel Ee{r,r') is given by the Fourier transform of the S'-matrix, 


Ee{r, r') 


^ f hl~^\kr)[l-Se{k)]h^,'^\kr')dk 

zn J 


Nb 

+ J2Anih^^\bnr)h^^\bnr') . 

n=l 


( 9 ) 


In the above equation, h^^\z) is the Riccati-Hankel function and A^e are the 
asymptotic normalization constants for the bound states E^ = —fi^/{2p)b\. 
The solution of Eq. (8) can then provide us the underlying local E-independent 
interaction V{r). 


f7(r) 


dKt{r,r) 

dr 


( 10 ) 


More details and practical aspects of the method can be found in Ref. [36]. 
We only stress here that the Marchenko equation is fully equivalent to the 
Schrodinger equation. However, while in the latter equation one uses the 
knowledge of the potential from which the physical information of the un¬ 
derlying system can be extracted via the wave-function, in the Marchenko 
equation the information at hand is the scattering and bound state data em¬ 
bedded in the S'-matrix. 

An efficient way to study systems with A > 3 is to employ the IDEA method 
to the A-body. The method is based on Hyperspherical Harmonics and the 
A-body wave function, due to pairwise acting forces, can be written as a sum 
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of Faddeev-type amplitudes for the pairs (f, j), 


'J'(x) = 51 = hf[L^](x) 5] F{rij,r), 


( 11 ) 


where is a harmonic polynomial of minimal degree Lm and r is the 

hyper-radius r = {2/A)J2i<j<A^ij- Using (11) one obtains the Faddeev-type 
equation for the amplitude Fijij^r) 

(r-B)%„,(x)f(ry,r) = -V^(r«)%„,(x) ^ F(r*,,r). (12) 

k<l<A 


This equation may be modihed by introducing in both sides the average of the 
potential over the unit hypersphere, i.e the so-called hyper-central potential 
Vo(r) [31,32] to obtain an integro-differential equation for the (modihed) Fad- 
deev amplitude (j){rij,r). Furthermore, assuming that the radial motion and 
the orbital motion are nearly decoupled, we may write 

0(rij,r) = P{z,r)/r^^~^^F ~ P^{z,r)ux{r) 

with 2 ; = 2r‘fj/r^ — l and one obtains the adiabatic approximation in which one 
extracts the potential surfaces U\{r) that provide us with the binding energies 
and the scattering states of the system. 

In contrast to the IDEA which is based on the Faddeev decomposition of the 
wave fnnction, in the AMD method [33,34] one assnmes that the antisym¬ 
metrized wave fnnction for the A-nncleons is given by a Slater determinant 

|t(Z)) = (1 ± F)(l/Vl!) det l^z.(rj)] (13) 


constrncted from the the i-th nucleon wave fnnction 


<^z,(d) 


(D#z,) \Xi) 


/ 2 z/^ 

3/4 

/ Zj 

2 lr. 2 ' 

V TT y 

1 exp 


1 + 2 ^* 


(14) 


In the above eqnations P is the parity projection operator, Z is a set of 
(complex) parameters parameterizing the nncleon, Z = (Zi, Z 2 , Z 3 ,..., Z^), 
z/ is a width parameter for the wave packet and \xi) = lo’jTj) is the spin- 
isospin fnnction for the nncleon i. To obtain the solution one employs the 
time dependent variational principle 


'/ {'l>(Z)|[tf4-ffl|1>(Z)) 

5 / dt - ^ 

*1 


($(Z)|$(Z)) 


= 0 


(15) 
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p (fm 

^He ^He ^<^0 

0 

8.67 30.7 1630 

0.003 

3.84 18.1 1280 

0.007 

0.75 8.56 979.6 

0.009 

— 5.70 874.8 

0.017 

— — 628.0 

0.034 

— ~ 507.8 


Table 1 

Binding energies (in MeV) in the eaa for ^He, ^He and systems for different 
densities p (in fm“^) 

where H is the total Hamiltonian of the system [33,34], 


One of the main advantages of the AMD formalism is the possibility of tracing 
the momentum and position of the individual wave packet. For this one dehnes 
the harmonic oscillator annihilation operator 


a = 




with 




( 16 ) 


where Zj is a complex vector with its real and imaginary component are related 
to the average position and moment of the Ath wave packet, 


Zi — + 



( 17 ) 


where 


D, 


(0Zi|0Zi) 


and Kj 


(0zJP#Zi) 

(0zj0Zi) 


(18) 


In our calculations we employ the general form of the nucleon-nucleon poten¬ 
tial 

V{r) = {W + BP, - HP, - MP,Pr)V{r) (19) 


where the bare potential V{r) is the Volkov [37] interaction which is widely 
used in model nuclear structure calculations. We first employ the Wigner force 
only to study the influence of correlations on the bound states of ^He, ^He and 
of for a various densities using the IDEA method. The binding energies 
obtained in the extreme adiabatic approximation (eaa) are given in Table 1. 
The ^He is the hrst to dissolve in the medium while the still is bound 


6 







Fig. 1. The two-body potential and the potential surfaces (in the extreme adiabatic 
approximation) for A = 2>, A = A and A = 16 for different densities p. 

at the quite high density of p = 0.034fm“^. Note, however, that already the 
isolated is strongly over-bound in the eaa. The disappearance of bound 
states is attributed to the weakening of the attraction in the interaction with 
increasing density p of the medium. 


This can be seen in Fig. 1 where the potentials surfaces (in the eaa) are shown 
together with the two-body interaction. The weakening of the force can be 
seen also in Fig. 2 where the Faddeev components P(r, z) are plotted for the 
“^He case. For p = 0.009 fm“^ the component is quite flat with the maximum 
shifted to ~ 5fm implying that the nucleons involved are already quite far 
from each other. 

.4 = 4 .4 = 4 

p = 0.009 




Fig. 2. The Faddeev components for the A = 4 system for p = Ofm ^ and for for 
p = 0.009 fm“^. The latter is just below the Mott density. 
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Fig. 3. Expectation values (r) for the center of the wave packets for protons (o) 
and neutrons (a) for the ®He, left figure, and ^Li, right figure, with a Wigner force 
only. The three configurations (up to down) correspond to p =0fm“^, 0.009 fm“^ 
and 0.017 fm“^ respectively. 
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Fig. 4. Expectation values (r) for the center of the wave packets for protons (o) and 
neutrons (a) for the ®He with Majorana force included. The three configurations 
(up to down) correspond to p =0fm“^, 0.009 fm“^ and 0.017fm“^ respectively. 

As already mentioned, the main advantage of the AMD method is the possibil¬ 
ity of tracing the movement of an individual proton or neutron by calculating 
the expectation value (r) for the corresponding wave packet. In Fig. 3 we show 
the results for the ®He and ^Li clusters with the Wigner force only. In the ^He 
case the two neutrons are squeezed out of the cluster leaving the a-particle 
behind which, however, also is enlarged and ready to fall apart. The ^Li is 
at hrst quite compact but two neutrons and a proton start to escape with 
increasing p and already at p = 0.017 fm“^ the 7 particles are not bound. It is 
interesting to note that although the a particle appears to be in the middle, 
the other particles always arranged themselves in an a particle formation. 

The effects of the inclusion of the Majorana force in the calculation is shown 
in Fig. 4. The differences generated by the inclusion of this force is striking. 
The a particle acts as an attractive center for the two neutrons forming a 
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halo nucleus which are escaping with increasing p. In the case where only 
a Wigner force is considered, the two neutrons are escaping in either side of 
the a particle. However, the inclusion of spin-isospin component in the nuclear 
force generates two sub-clusters consisting of an a cluster and the two neutrons 
which tend to split hrst prior an overall disintegration occurs. These differences 
are best seen on the two-dimensional Fig. 5 obtained by integrating over the 
2 :-axis. 

Our conclusions can be summarized as follows, i) The dynamical behavior 
of a nucleonic cluster in a medium can be reliably studied using a combi¬ 
nation of methods namely, a cluster mean held approximation of the Dyson 
Green function expansion, an inverse scattering method for transforming the 
£^-dependent potential to an i?-independent one and few-body methods, ii) 
The construction of the i?-independent interactions paves the way to study 
the behaviour of light nuclei (A~3-20) in a medium using conhguration space 
formalisms, iii) The use of the AMD method allows us to trace the individual 
movement of the nucleonic wave packet and the formation of sub-clusters when 
the density increases. This feature is particularly suited to explore correlations 
as has been done, e.g., for light clusters in a recent experiment [8], and also 
get insight into a possible evaporation mechanism of light clusters, iv) The 
use of Wigner forces alone seems to be insufficient for a proper description 
of the clusters while the inclusion of the Majorana force generated dramatic 
changes in the clustering picture of the particles. However, in order to obtain a 
complete picture on the behavior of quasi-nuclei in a changing nuclear matter 
environment, more realistic forces that include also the Bartlett and Heisen¬ 
berg components as well as an LS-force should be employed. Although the 
effects of the former two components are expected to be small, the role played 
by the latter force in these extremely loose bound systems could be quite 
important and therefore its inclusion in studies concerning disintegration of 
clusters is warranted. 
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Fig. 5. Two-dimensional density plots for the ®He system with a Wigner force (left 
panel) and with Majorana force included (right panel) for p = Ofm“^, 0.007fm“^ 
and 0.017fm“^ (top to bottom). The figure scales are in fm. In the left panel the 
outermost line corresponds to the value of 0.01 fm“^ and gradually increases by 0.5 
fm“^ until it reaches the value of 5 fm“^ for the inermost line. In the right panel 
the correponding numbers are 0.01 fm“^ until 2 fm“^ in steps of .2 fm“^. 
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